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ABSTRACT 



The correlation between the two r-polarisations is predicted to be opposite (H — / — h) for 
the H/A — > t + t~ signal, while it corresponds to the same sign combinations for the Drell-Yan 

(+ + / ) and tt ( ) backgrounds. We show that this correlation can serve as a distinctive 

test to confirm the presence of the MSSM Higgs bosons H and A in their hadronic r + r~ decay 
channels at the LHC. 
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The Minimal Supersymmetric Standard Model (MSSM) contains two Higgs doublets [1]. The 
ratio of their vacuum expectation values is denoted by tan/3, which is required to be larger than 
2 by the LEP data [2, 3]. The two complex scalar doublets correspond to eight independent 
fields, three of which are absorbed by the W ± and Z bosons. Thus we are left with five physical 
states: two neutral scalars h and H, a neutral pseudoscalar A, along with a pair of charged 
Higgs bosons H ± . For large Ma values (>> Mz), the two heavy neutral Higgs bosons, H and A, 
have similar masses and couplings. On the one hand, the H/A boson couplings to the top quark 
are suppressed by a factor of cot/9 relative to that of the Standard Model (SM) Higgs boson, 
while they both decouple from the W ± and Z bosons. On the other hand, their couplings to 
the bottom quark are enhanced by a factor of tan/5 relative to that of the SM Higgs boson. 
Thus one expects a large H/A production signal at the Large Hadron Collider (LHC) over the 
large tan/5 region via the bb — > H/A 'fusion' mechanism. Moreover the H/A couplings to the r 
lepton are also enhanced by a factor of tan/3 over this region. 

Therefore the H/A — ► t + t~ decay channel provides the most promising signature for the 
two heavier neutral Higgs bosons of the MSSM at the LHC. Simulation studies have shown 
this decay channel to be viable over a large part of the parameter space of the model [2, 3]. 
Moreover, a systematic study of this process in both the leptonic and hadronic decay channels 
of the r's [4, 5] shows that the latter has two advantages, i.e., a larger branching ratio (BR) 
and a larger visible fraction of the r-energy/momentum. This implies that one can reconstruct 
the r + r _ invariant mass most precisely via 

t + t- -> h + h~vv (1) 

(hereafter, h+ represents a generic charged hadron) [4]. This channel is particularly useful for 
relatively large values of the two MSSM parameters defining the Higgs sector of the MSSM at 
tree level, i.e., M4 ^ 200 GeV and tan/5 ^ 15. This area is precisely where the two Higgs 
states H and A are degenerate in their masses and couplings. 

In this note, we shall analyse a distinctive feature of the H/A — > t + t~ signal that has 
not received much attention so far, i.e., the correlation between the r-polarisations. It simply 
follows from angular momentum conservation that the t + t~ pairs emerging from the H/A 
decays have opposite polarisations, i.e., 

H/A - t+t-X, P T+ (P T -) = ±l(Tl). (2) 

In contrast, the main irreducible backgrounds from Drell-Yan (DY) and top-antitop processes 
predict the same polarisation for the t + t~ pair, i.e., 

7 */ Z ^ T + T - X , P T+ (Pr-) = ±1(±1), (3) 

tt -> bbW + W~ -> t+t- X, P T+ (P T -) = -l(-l). (4) 

The distinctive r-polarisations for the above signal and background processes were first dis- 
cussed in Ref. [6], along with the analogous distinction between the charged Higgs signal and 
the W+ background, i.e., 

H± -> t±v, P T ± = +1; W± -> t±v, P T ± = -1. (5) 
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The possibility of measuring the predicted r-polarisations via their hadronic decays was also 
discussed there. Meanwhile, quantitative Monte Carlo (MC) studies of the charged Higgs boson 
signal at the LHC have shown the polarisation distinction in (5) to be a very powerful tool in 
separating the if 1 * 1 and W ± contributions in the hadronic r-decay channels [7, 8]. To the best of 
our knowledge, however, there has been no such quantitative study for the neutral Higgs signal 
in (2). Of course, one can see from Eqs. (2)-(4) that the polarisation distinction for the neutral 
Higgs signal is more subtle than the charged Higgs case, Eq. (5). In particular, the average 
r-polarisation for the H/A — > r + r~ signal is 0. The same is also true for the DY background 
in (3), since the vector coupling of 7 and the dominantly axial coupling of Z imply that the 

polarisation combinations ++ and occur with almost equal probability. Nonetheless, we 

shall see below that the correlation between the r-polarisations can be effectively used to 
distinguish the signal from the background in the hadronic r + r~ channel. 

Our results are based on a parton level MC simulation of the H/A signal in (2) as well as 
the 7*, Z and ti background processes in (3) and (4). The dominant signal process in the large 
tan/3 (<; 15) region of interest can be emulated via associate production with 6-quarks [9]: 

gg,qq^bbH/A. (6) 

The cross sections for process (6) have been calculated by several groups [10]. While these 
references do not give compact formulae, the latter are very similar to those for gg, qq — > tbH~ , 
described in detail in Ref. [11], which have been used again here, after rearranging the Higgs-to- 
fermion couplings appropriately. The formulae for the H/A decay rates are taken from Ref. [12] 
(see also [13]). The DY and top-antitop background processes are computed using MadGraph 
[14]. The phase space integrations have been done numerically with the help of VEGAS [15]. 
We have used the parton distribution functions MRS99(COR01) [16], evaluated at the scale 
fj, = Mh(Ma) for the signal and /x = V§/2(2m t ) for the 7*, Z{tt) backgrounds. (The results for 
these two processes are however insensitive to the precise value of these scales.) The 'running' 
expression of the 6-quark mass was used for the Yukawa coupling of the signal process. 

To emulate detector resolution, we have applied Gaussian smearing on the transverse mo- 
mentum of each parton jet (including the r-jets) with 

(<j( Pt )/pt) 2 = (O.e/^) 2 + (0.04) 2 (7) 

in GeV units. The missing transverse momentum, p™ lss , is reconstructed from the vector sum 
of the jet pt's after resolution smearing. Following Ref. [4], we have imposed the following 
selection cuts on the two r-jets and the accompanying missing transverse momentum: 

Ai?(jet - jet) > 0.4, p T (r-jets) > 60 GeV, |r/(r-jets)| < 2.5, 
A0(r-jets) < 175°, p^ iss > 40 GeV. (8) 

As in Ref. [4], we shall concentrate on the 1-prong hadronic decay channels of r's, which are 
best suited for distinguishing r-jets from the QCD background. With the above cuts, the QCD 
background has been estimated to be relatively small compared to the DY and top-antitop 
ones [4]. While the DY background could be suppressed by 6-tagging, it would also reduce the 
signal size significantly. Therefore, we do not require any 6-tagging, just like Ref. [4]. 
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The 1-prong hadronic decay channel of r's accounts for 80% of its hadronic decay width 
and 50% of the total width. The main contributors to the 1-prong hadronic decay are 



r ±-> 7^(12%), p ± z/(26%), afv(8%), (9) 

where the BRs for n and p include the small K and K* contributions, respectively, which have 
identical polarisation effects [17]. The centre-of-mass (CM) angular distribution of r-decays 
into 7r or a vector meson v(— p, a±) is simply given in terms of its polarisation as 

1 " ! hi + P rC0S 9), (10) 



T n dcos9 2 
1 dT vLjT _ \m\,m\ 

where L, T denote the longitudinal and transverse polarisation states of the vector meson. The 
fraction x of r-momentum in the laboratory frame carried by its decay hadron is related to the 
polar angle 6 via 

1 "IT? 

x= _(l + O os0) + ^ (l-cosfl), (12) 

in the collinear approximation. Thus the visible momentum of r, which defines the momentum 
of the r-jet, is 

Pr-jet = Phadron = Xp T . (13) 

We see from Eqs. (10)-(13) that P T = +1 gives a harder r-jet than P T = —1. 

We can also see from Eqs. (10)— (13) that the hard r-jet is dominated by the n and the 
longitudinal vector meson (pl,o>il) contributions for P T = +1, while it is dominated by the 
transverse p and a\ (pr, air) contributions for P T = —1. The two sets can be distinguished 
by exploiting the fact that the transverse p and a\ decays favour even sharing of momentum 
between the decay pions, while the longitudinal p and a\ decays favour uneven sharing, where 
the charged pion carries very little or most of the meson momentum. One can measure the 
fraction of the visible r-jet momentum, carried by the charged pion, 

R = Ptt± /Pr-jet, (14) 

by combining the charged prong momentum measurement in the tracker with the calorimetric 
energy deposit of the r-jet. Then the hard r-jet is predicted to show peaks at R ~ and 1 
for P T = +1 and at R ~ 0.4 for P T = —1 [6, 7]. Note that the r^ — > ix^v contribution adds 
substantially to the R ~ 1 peak for P T = +1, while it is suppressed for P T — — 1. One can 
easily derive the above results quantitatively for pi,r decays; but one has to assume a dynamical 
model for the a\ decay for a quantitative result. 

We shall adopt the model of Ref . [18] , based on a conserved axial- vector current approxi- 
mation, which describes the a\ — > 3n decay data very well. A detailed account of the p and 
ai decay formalisms including finite width effects can be found in Refs. [6, 7]. A simple FOR- 
TRAN code for 1-prong hadronic decays of polarised r-leptons based on these formalisms can 
be obtained from one of the authors. 
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Before discussing the results, let us briefly describe the reconstruction of the t + t" invariant 
mass. The A0 cut of Eq. (8) ensures that the transverse momenta of the two r-jets, pn and 
p T 2: are not in the back-to-back configuration. In this case, one can resolve the p™ lss along their 
directions and add it to pn and pT2- Scaling up the respective sums by the ratios pi/pri and 
P2/PT2 gives the reconstructed momenta of the r + r~ pair. The resulting t + t~ invariant mass 
represents the H/A mass for the signal process (2) and the 7*/^ mass for the DY background 
(3); but it does not represent any physical quantity for the tt background (4), since the latter 
contains additional sources of missing momentum. 

The production cross sections of the H and A bosons are practically identical in the large 
tan/3 region of our interest. As an illustration, we have calculated the signal cross sections for 

tan/3 =15 and M4 = Mh = 200(300) GeV, (15) 

which is at the edge of the discovery limit for this channel [4, 8]. 

Fig. 1 shows the signal and background cross sections against the reconstructed r + r~ in- 
variant mass. We see that a M TT cut of 150-300(200-400) GeV will suppress the DY and 
top-antitop backgrounds without any appreciable loss of the 200(300) GeV H/A signal. There- 
fore, we shall impose these invariant mass cuts in the reminder of our analysis. Let us first 
consider the Mh,a = 200 GeV case. The signal and background cross sections for this case are 
shown on the top row of Tab. 1. They are seen to be of similar size. With an expected lumi- 
nosity of 300 fb" 1 at the LHC, they correspond to about 1000 signal events over a background 
of roughly similar size. This constitutes a comfortably large S/ \f~B ratio. 



M HA (GeV) 


bbA 


bbH 


7 *,Z 


tt 


200 
300 


1.58 
1.88 


1.57 
1.86 


3.03 
4.14 


0.73 
0.60 


tan/3 = 15 



Table 1: Cross sections in femtobarns for Higgs signals and backgrounds in the r + r X channel, 
after the selection cuts of Eq. (8). In addition, we have imposed a M TT cut of 150-300 GeV for 
M HA = 200 GeV and 200-400 GeV for M HA = 300 GeV. We have also imposed a stronger p T 
cut of 100 GeV on the harder r-jet in the latter case. 

Fig. 2 shows the three-dimensional (3-D) plots of the 200 GeV H/A signals along with the 
DY and tt backgrounds against R\,R2, the fractional r-jet momenta carried by the charged 
prongs. The subscript 1(2) refers to the harder (softer) of the two r-jets. There are several 
points worth noting in this figure. The r-identification will require a minimum hardness for 
the charged prongs, presumably R\ and R 2 > 0.2. Thus the R± j2 < 0.2 contributions from p L 
and an, will not be relevant for the practical analysis. However, the peaks at R±^ — 0.4 and 
R\,2 — 1, coming from the pt, o>it and Pl,^ respectively, can be used to distinguish between 
the P T = — 1 and +1 contributions. In particular, the tt background, corresponding to the 
polarisation combination , is crowded around R\^ — 0.4. In contrast, the signal is peaked 
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at Ri or R 2 > 0.8, as expected for the polarisation combination H — . The same is seen to be 

true for the DY background, because the px > 60 GeV cut on the r-jets has suppressed the 

contribution with respect to the ++ combination. Thus requiring at least one of the two r-jets 
to contain a very hard charged prong, carrying > 80% of its visible momentum, will suppress 
the ti background effectively without any significant loss to the signal or the DY background. 
Moreover, we expect this requirement to suppress the QCD background to an even greater 
extent than the ti background, because its R distributions are softer than those of the latter 
[4]. Therefore, we strongly advocate an asymmetric cut, requiring at least one of the Ri and 
R 2 to be > 0.8, for effective suppression of the tt and QCD backgrounds (including W + jets). 
Then the signal can be distinguished from the remaining DY background by looking at the 
asymmetric region: Ri > 0.8 and R 2 = 0.2 — 0.8 and vice versa. The signal(DY background) 
has a peak(dip) in this region as expected from the polarisation correlation H — (++)• We have 
estimated 55% of the signal events to come from this asymmetric region as against 43% for the 
DY background. Thus excluding the symmetric region will improve the S/B ratio, but not the 
S/ \J~B. On the other hand, the distribution of events between the symmetric and asymmetric 
regions can be used to distinguish the signal from the DY background. From an event sample 
of ~ 1000 one expects 550 to populate the above asymmetric region for the signal (H — ) as 
against 430 for the background (++), which constitutes a ^ 5a excess. 

For higher H/A masses it may be more advantageous to use asymmetric pt cuts on the 
two r-jets. Fig. 3 shows the 300 GeV H/A signal along with the corresponding DY and ti 
backgrounds against the pr of the harder r-jet. It shows that increasing the pn cut to 100 GeV 
(while keeping the pti cut at 60 GeV) will suppress the tt background to a much higher degree 

than the signal, as expected from their polarisation combinations, and — h, respectively. 

The DY background is also suppressed to a somewhat larger extent than the signal, though the 
difference in rather modest (7-8%) in this case. The pT2 > 60 GeV cut is already suppressing 

the contribution with respect to the +- 1- combination even for M TT ~ 300 GeV. We expect 

this difference to increase for higher H/A masses. We also expect the higher p T \ cut to be very 
effective in suppressing the QCD background (again, including W^+jets). Tab. 1 shows the 
300 GeV H/A signal cross sections along with the DY and tt backgrounds for this asymmetric 
cut (i.e., pti > 100 GeV, pr2 > 60 GeV). Their sizes are similar to the 200 GeV signal and 
background cross sections with symmetric pr cuts. We expect the asymmetric pr cuts to be 
advantageous in probing the region M A / H > 300 GeV. 

Fig. 4 shows the 3-D plots of the 300 GeV H/A signal along with the DY and ti backgrounds 
against R\ and R 2 for the asymmetric px cuts mentioned above. Most of the observations 
made for Fig. 2 apply to this case as well. In particular, one sees a peak(dip) for the signal(DY 
background) at Ri > 0.8 and R 2 = 0.2 — 0.8, as expected from the polarisation correlation 
H — (++)• Note, however, that there is no clear peak for the signal in the complementary 
region, R 2 > 0.8 and i? x = 0.2 — 0.8. This is because in this case the p T 2 > 60 GeV cut 
is not hard enough to suppress the pion peak (at R 2 — 1) for P T2 = —1. Consequently, the 
accompanying T\ has a significant component of P Tl = +1. Thus for higher H/A masses one 
should look for the signal peak in the Ri > 0.8 and R 2 = 0.2 — 0.8 region only. Nonetheless 
the difference is very significant in this case. About 30% of the 300 GeV H/A signal events 
populate this asymmetric region against 20% of the DY background events. This constitutes a 
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7(T excess for an event sample of ~ 1000, which can be used as a confirmatory test of the H/A 
signal. 

Finally, notice that the percentages of the signal and background events populating the 
asymmetric region are independent of the normalisation uncertainties of the various processes, 
since they depend entirely on the polarisations of the latter. The former mainly stem from 
higher order corrections due to QCD. While these are relatively well under control for the two 
backgrounds considered here, they can be very large for the signal (see [19] for an up-to-date 
discussion). 

In summary, we have found that the correlation between the r-polarisations are likely to play 
a very useful role in distinguishing the H/A signal from the backgrounds in the hadronic t + t~ 
channel at LHC. In view of the importance of this channel in the MSSM Higgs search programme 
at the LHC, it is imperative that this analysis is followed up with a more rigorous MC simulation, 
incorporating jet hadronisation and detector acceptance effects. We have currently undertaken 
such a simulation study [20] using the HERWIG event generator [21, 22] interfaced with the 
TAUOLA package [23, 24] for polarised r-decays. 

We gratefully acknowledge discussions with P. Aurenche, D. Denegri., R. Kinnunen and A. 
Nikitenko. 
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Figure 1: Singly differential distributions in the reconstructed t + t~ invariant mass in 1- 
prong hadronic decays, for the Higgs signals (dashed and dotted lines for H and A, respec- 
tively, coinciding visually) and backgrounds (solid and dot-dashed lines for DY and top- 
antitop, respectively) in the t + t~X channel, after the following selection cuts (at parton 
level): Pt{t — jetmi) > 60[60] GeV (1[2] refers to the mostfleast] energetic), \t)(t — jet) | < 2.5, 
A$(t - jets) < 175°, Ai?(jet-jet) > 0.4 and pf iss > 40 GeV. Normalisations are to the total 
cross sections at y/s = 14 TeV (for Ma = M H = 200(300) GeV and tan/3 = 15, in the case of 
the signal). Bins are 10 GeV wide. 
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Figure 2: Doubly differential distributions in the energy fractions Ri = pj jp\ ~ jct , where p^ 
is the charged pion momentum and p\ ~ Jct that of the visible r— jet, i.e, the momentum carried 
away by the mesons 7r ± ,p ± and af in 1-prong decays, for Higgs signals and backgrounds in 
the t+t~X channel (with pi > p 2 ), after the following selection cuts (at parton level): prij — 
jet 1[2] ) > 60[60] GeV (1[2] refers to the most [least] energetic), ^(r — jet)| < 2.5, A$(r— jets) < 
175°, Ai?(jet-jet) > 0.4, pf ss > 40 GeV and 150 GeV < M TT < 300 GeV. Normalisations are 
to unity at ^/s = 14 TeV (for M A = M H = 200 GeV in the case of the signal). Bins are 0.2 
units wide. 
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Figure 3: Singly differential distributions in the transverse momentum of the most energetic r- 
jet in 1-prong decays, for Higgs signals and backgrounds in the t + t~X channel, after the 
following selection cuts (at parton level): p T (r — jet^]) > 60[60] GeV (1[2] refers to the 
most [least] energetic), Mr -jet) | < 2.5, A$(r-jets) < 175°, Ai?(jet-jet) > 0.4 and^ iss > 40 
GeV. Normalisations are to unity at ^Js = 14 TeV (for Ma = Mh = 300 GeV in the case of 
the signal). Bins are 4 GeV wide. 
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d 2 a/dR l /dR 2 /a 



bbA -> t+t-X (M A = 300 GeV) bbH -> t+t~X (M h = 300 GeV) 





Figure 4: Doubly differential distributions in the energy fractions Ri = pj /Pi~ iC , where 
pj is the charged pion momentum and p[ _Jct that of the visible r— jet, i.e, the momentum 
carried away by the mesons ti ± ) and af in 1-prong decays, for Higgs signals and backgrounds 
in the t + t~X channel (with p\ > p 2 ), after the following selection cuts (at parton level): 
Pt(t — jet 1 [ 2 ]) > 100[60] GeV (1[2] refers to the mostfleast] energetic), \t](t — jet) | < 2.5, 
A$(r - jets) < 175°, Ai?(jet-jet) > 0.4, pf ss > 40 GeV and 200 GeV < M TT < 400 GeV. 
Normalisations are to unity at y/s = 14 TeV (for Ma = M H = 300 GeV in the case of the 
signal). Bins are 0.2 units wide. 
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